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Shadows of multi-black holes have structures distinct from the mere superposition of the shadow 
of a single black hole: the eyebrow-like structures outside the main shadows and the deformation 
of the shadows. We present analytic estimates of these structures using the static multi-black hole 
solution (Majumdar-Papapetrou solution). We show that the width of the eyebrow is related with 
the distance between the black holes and that the shadows are deformed into ellipses due to the 
presence of the second black holes. These results are helpful to understand qualitatively the features 
of the shadows of colliding black holes. We also present the shadows of colliding/coalescing black 
holes in the Kastor-Traschen solution. 

PACS numbers: 97.60.Lf ; 04.70.-s 

I. INTRODUCTION 

One of the hottest topics in galaxy formation is the coevolution of super massive black holes (SMBHs) with spheroid 
components (bulges) of galaxies. It becomes more and more clear that most of galaxies and AGNs have at least one 
SMBH and there is a strong correlation between the SMBH mass and the bulge mass of host galaxies 

PH. Although a 

detailed mechanism of the coevolution is not yet understood, it is almost certain galaxy mergers play an essential role, 
since it is known that bulges or spheroid components are formed by merger of galaxies in the hierarchical clustering 
scenario of structure formation. Hence it may be also natural to consider formation of SMBHs is due to mergers of 
smaller black holes. 

An observation clue of existence of binary black holes is recently obtained from detailed study of Kepler motion of 
a radio emission component in the radio galaxy 3C 66B by using a technique of phase-referencing very-long-baseline 
interferometry (VLBI) In particular, a newly found periodic flux variation suggests that this binary system will 
coalesce in 500 years [5(. Perhaps, wc may conclude coalescence of binary black holes often takes place in the Universe. 
However a direct evidence of black hole merger is still missing. One of the possibilities to see the merger process is to 
observe shadows of black holes shone by the radiation from the accretion disc or star lights behind the black holes [f| . 
Since two event horizons merge into one event horizon, we expect that the shadows must show very peculiar time 
evolution. Observing these shadows, therefore, should be compelling evidence of a coalescing black holes as well as 
provides an intriguing probe of general relativity with very strong gravity field. 

As a first step toward the study of a realistic black hole binary we have recently calculated the shadows of the 
Kastor-Traschen (KT) 0] cosmological multi-black hole solutions [gj. We have found that the shadows arc deformed 
in the direction of the collision and that in addition to the shadow of each black hole, eyebrow-like structures appear 
as the black holes come close to each other. In this paper, we attempt to understand these structures analytically 
using the Majumdar-Papapetrou (MP) solution Q, the static multi-black hole solution with the maximal charge to 
which the KT solution is reduced when the cosmological constant is zero. 

The paper is organized as follows. In Sec. II, we present several analytic calculations of null geodesies in the 
extreme Rcissncr-Nordstrom and MP solutions in order to provide analytic estimate of the eyebrow-like structures of 
the shadows as well as the deformation of the shadows in the MP solution. In Sec. Ill, extending our previous results 
Q, we present numerical results of the shadows of the colliding/coalescing black hole binary in the KT solution by 
changing the viewing angle of the distant observer and compare them with the shadows of the MP solution. Sec. IV 
summarizes the results. 
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II. SHADOWS OF STATIC TWO BLACK HOLES: MAJUMDAR-PAPAPETROU SOLUTION 
A. Geodesies in Extreme Reissner-Nordstrom solution 

The Majumdar-Papapetrou (MP) solution Q is given by 

ds 2 = -£l- 2 dt 2 + il 2 (dx 2 + dy 2 + dz 2 ), (1) 

r> = l + V— , n = ^(x - Xi y + (y- yi ) 2 + (z- ztf, 
^ n 

% 

where m, is the mass of i-th black hole located at r = (we use the geometrical units, G = c = 1). The MP solution 
in the case of a single black hole is reduced to the extremal Reissner-Nordstrom (RN) solution. 

First of all, we derive the relation between the deflection ang le and the impact parameter close to the unstable 
circular orbit for the extreme RN solution following Luminet |lOl ITT1 ] , in order to understand the shape of black hole 
shadows in MP spacetime. 

The metric of the extreme RN solution is given by 

ds 2 = - (l - 2 dt 2 + (l - 2 dR 2 + R 2 (d9 2 + sin 2 6d0 2 ) , (2) 

and the horizon is at R = M. The radial coordinate R is related with the isotropic coordinate r in Eq. (JTJ) as 
R = r + AI. From the spherical symmetry, we may restrict ourselves to the equatorial plane without loss of generality. 

In terms of two conserved quantities, i.e. the energy E = ^1 — —j i and the angular momentum L = R 2 4>, the null 
geodesies satisfy the "energy equation" 

1 -R 2 + V{R)= 1 -E 2 , (3) 
L 2 ( Ai^ 2 



where t = dt/dX with A being the afhnc parameter. V(R) is shown in Fig. [T] The null geodesies with E 2 /2 being 
larger than L 2 /32M 2 , which is the local maximum of V(R) at R = 2M, fall into the black hole. In terms of the 
impact parameter b = L/E, the light rays with b < 4M fall into the black hole. Thus, the shadow of the extreme RN 
black hole is the disk of radius AM. 

Let us investigate the behaviors of the null geodesies approaching the unstable circular orbit at R = 2M. In terms 
of u = 1/R, Eq. ([3]) is written as 

f^j 2 = -u 2 (l-Muf + ^f(u). (4) 

For b > AM, the quartic equation f(u) = has four distinct roots, three of which are positive. Arranging the ordering 
of the roots as u\ < < U2 < < u^, l/u2 corresponds to the periastron distance. Then the solution of Eq. (j4j can 
be written as 

du 




b + AM 



K{k)-F[^k) , 



where the origin of <$ is chosen at the periastron passage when u = u 2 and k is given by 



k 2 (6) 
k ~ b + AM ' {b) 



F((p, k) is the elliptic integral of the first kind defined by 

F(<p,k)= / - (7) 
Jo VI - k 2 sin 2 6 



3 




FIG. 1: The effective potential V(R) for the extreme RN black hole. 



and K(k) = F(ir/2, k) is the complete elliptic integral of the first kind. 

The solution can be used to obtain the asymptotic behavior of for l/it2 — > 2M. If we write \ju% = 2M(1 + 8) 
with 8 <C 1, then from /(U2) = in Eq. 6 = 4A/(1 + J 2 ) and hence k 2 = 1 — 4p Using the asymptotic relation 
*T(fc) -> § In (16/(1 - A: 2 )) for fc -> 1, we obtain 

1 2 5 

4>oo = -f= In — f= ■ (8) 

^ V2 ^/2 + l) 2 <5 2 1 ; 

Therefore, the light rays with the impact parameter b = 4M(1 + <5 2 ) are deflected by the black hole by the angle 
= 20oc — 7r with the relation 

b = UI + — 2?M e~^e~75 e . (9) 
(V2 + 1) 2 

For those geodesies that go round the black hole n-times, the deflection angle can be written as + 2wr. Accordingly, 
we have 

2 7 M 



4A/+ ZM 2 e-^e-^ {e+2n *\ (10) 
(V2 + 1) 



B. Shadows of MP solution 



Next, we study the shadows of two black holes in the MP spacetime. We consider a two equal mass black hole 
system in the MP spacetime Eq. ([T]). The black holes are located at z = ±£/2 with mass m = M/2. Because of 
the axial symmetry, there exists a conserved quantity L z (corresponding to the angular momentum integral) for the 
geodesies as well as the energy integral which are defined in the spherical coordinate by 

E = n~ 2 i, L z = fl 2 r 2 sin 2 6<j>. (11) 

We calculate the null geodesies from the observer located far from the black holes. The geodesies that fall into the 
black hole form the shadow of the black holes. In Fig. [2J the shadows of two black holes in the MP spacetime are 
shown. The observer is located at the equatorial plane (6 = tt/2). Here we have defined the celestial coordinate of an 
observer, (a,/3), as 

r pW r p(e) 
a = ^ — pWT' & = ^ ~pWT> ( 12 ) 



4 



-12 -8 -4 4 8 12 -8 -4 4 8 




-8 -4 4 8 -8 -4 4 8 



J3/M 

FIG. 2: Shadows of two black holes in the MP spacetime for several distances t between the black holes. 

where are the momenta in the local inertial frame. 

We find that there appears additional eyebrow-like structure in the outer region of the main shadow. The eyebrow 
grows as the distance between the black holes is closer. Although not discernible in the figure, in fact there appears 
"the fine structure" of the eyebrows: infinitely many thinner eyebrows at the outer region of these eyebrows as well 
as at the inner region of the main shadow. 

Moreover, the distance between the main shadows is slightly larger than the real distance between the black holes, 
and the shapes of the shadows are suppressed in the (3 direction and are slightly elongated in the a direction. Note 
that the radius of the shadow of a single black hole is Ami. Take £ = 2M case as an example. If the other black hole 
is absent, the shadow should be located at /3/M = ±1 with the radius Am = 2M. However, because of the presence of 
the second black hole, the inner edge of the shadow is shaved and the outer edge is stretched toward larger /3 instead. 
In the following we study how these structures, different from the superposition of the shadow of a single black hole, 
appears. 

C. Eyebrow in Black Hole Shadow 

First, we try to explain how the "eyebrows" are formed. We restrict ourselves to the a — plane. The schematic 
picture of null geodesies for the silhouette of the eyebrows is shown in Fig. [3] The geodesies coming in with impact 
parameter b + and &-(< b + ) arc deflected by the left black hole. The geodesies with 6+ go into the unstable circular 
orbit of the right black hole from the outside, while inside for those with 6_. The geodesies with the impact parameters 
between 6_ and b+ form the eyebrow of the shadow. Hence, the width of the eyebrow Ab = b + — &_ is directly related 
with the distance between the hole £} 

Then we calculate the relation between the deflection angle and the impact parameter and the distance. We 
assume the distance is so large I 3> m that the effect of the second black hole on the spacetime may be treated as a 
perturbation to a single black hole. The mass of the black holes is assumed to be equal with m = M/2. From Fig. Q] 



Likewise, the fine structure of the eyebrows may be understood as those null geodesies that are captured by the one hole after going 
around the other hole several times. 




FIG. 3: Null geodesies for the "eyebrow" of black hole shadows. The geodesies coming in with the impact parameter b+ and 
&_(< 6+) are deflected by the left black hole. Then they go into the unstable circular orbit of the right black hole: the one with 
6+ from the outside, while inside with 6_ . The geodesies with the impact parameters between fe_ and 6+ form the eyebrow of 
the shadow. 





FIG. 4: The relation between b+, 0+ and I. 



FIG. 5: The relation between 6_, 0_ and I. 



and Fig. [5j wc find that the deflection angle Q± of the geodesies with the impact parameter b± is given by 

„ 7r b + — 4m „ 7r &_ + 4m 



Then from Eq. © 

A5 _ b+-b 
m m 



m (a/2 + 1) 

Hence we obtain to the first order in raj I 



e V2 I e — e ^ 



2 7 m 



Ab 



2 10 m 2 



(V2 + 1) V2£ 



-e 2V2 ~ i.i- 



2 + 1) y/2l 



M 2 



Ab 

m 



(13) 



(14) 



(15) 



We numerically compute the width of the eyebrow Ab for several I. The results are shown in Fig. [6l The least 
square fit of the numerical data for 10 2 < £/M < 10 4 gives 



In — = 0.96 - 0.99 In [ — I , 
M \ M 



(16) 



6 



Ab/M 



0.01 



0.001 



0.0001 




10000 



FIG. 6: The width of the eyebrow A6 as a function of the distance between the holes. The numerical results (dot: red) and its 
linear fit (solid: black). The dashed (blue) line is the analytic relation Eq. (|15[1 . 



or Ab/M = 2.6(£/M) °". Apart from the slight offset, the power index is in good agreement with the analytic 
estimate Eq. (fT5j) . 



D. Deformation of Black Hole Shadow 

Next, we study how the shadows are deformed and are shifted toward a larger separation. 

We assume that two black holes are located at z — and z — £. In the spherical coordinate, f2 is given by 



n = 1 + ^ + = l+ 'Jl+^O), r' = vV + £ 2 -2£rcos6. 



m 
r 



(17) 



We assume £ 3> m, m', in order to treat the effect of the black hole m! at z ~ £ as a perturbation to the black hole m 
at z = 0. 

In addition to the conserved quantities E and L z (Eq. we also dehne the quantity Q which corresponds to 

the Carter constant for a single black hole 



Q = P g ' + cot z 9L Z = L x + Ly, 
where Pq = Q 2 r 2 8. Introducing ( 2 = Q/E 2 , £ = L z /E, we have for the null geodesies 



(18) 



(19) 



Note that ( (or Q) is no longer conserved due to the presence of the black hole m! . In fact, from the null geodesic 
equation for 8, we have 



Ann oPo 



Introducing [i = cos 8, then in terms of the derivative with respect to r, Eq. (p~9j) and Eq. (|20| can be written as 



dr 

dr dr 



c 2 -(c 2 +e 2 y 2 
r 4 ^ 4 - {c 2 +ey 2 



(20) 

(21) 
(22) 
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The equations determining the unstable orbits are 

C 2 + £ 2 

K=l-^- = (23) 

dr ~ r 2 4 dr + r 4 Q 8 dr 1 j ~ U ' 1 j 

Eq. ([23]l. Eq. ((24]) and Eq. j22]| are combined to give 

c 2 +e 2 = r 2 fi 4 j (25) 

n + 2rfi r = 0. (26) 

Denoting the solution to the above equations as r c , r c describes the shape of the photon sphere deformed by the 
second black hole. 

We set to = 1 and to' = k and measure the length in units of m for notational simplicity. Then for a single black 
hole, the radius of the photon sphere r c is r c = 1 and the radius of the black hole shadow is 4. 
We expand C according to the order ofl/£ for fixed r, 

C = Co + Ci + C 2 + ■ • ■ • (27) 

The radius of the photon sphere r c is also expanded 

r c = l + e 1 +e 2 + .... (28) 

We consider the observer located at 6 ~ tt/2. Then rpi/£ <C 1 holds everywhere on the unstable orbit. So ip in Eq. 
(|17p can be expanded in terms of /.t 

V>(r, /i) = Vo(r) + Mr)n + \i>i(r)f + ^ipa Mm 3 + ..., (29) 

where 

k £r ( Ir \ " 

^ 0= V^TF' ^ 1=V VT^' ^=M2n-l)"(^+^J • (30) 

For the observer located at (r, 6, </>) = lim r _ i . 00 (r, 7r/2, 0), £ and C are related to the celestial coordinate a and (3 (Eq. 
(HU)) as £ = a, Co(r -> oo) = /3 |p. 
As explained in Appendix [A] it is shown that up to 0(l/£), ( satisfies the following equation 

§- = ( 31 ) 



(W/i order 



-2Vi, (32) 



From Eq. (|3T|) . Co satisfies 

where we have neglected 1/r as a small perturbation. Since Co ~~ ► P as r — ^ 00 > the solution is 

2k^ , , 

We need to determine r c to calculate Co- From Eq. (f2l))) . r c — 1 to the zero-th order and then Co = P + 2k. Hence 
using Eq. (|25[) . we find 

a 2 + (/3 + 2k) 2 = 16. (34) 

In the zero-th order, the shape of the black hole shadow is a circle with the radius of 4 which is the same as a single 
black hole. However the center is shifted toward j3 < direction by —2k. This explain why the distance between the 
main shadow is slightly larger than the real distance between the black holes. The reason is simple: the null geodesies 
with /3 > are attracted not only by the black hole at z = but also by the black hole at z = £ in the opposite 
direction, so they are deflected less, while for j3 < the null geodesies are attracted by both black holes in the same 
direction, so they are deflected more. Hence, this feature should not be limited to the MP spacetime, but should be 
present for uncharged black holes. 
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2. 1st order 



From Eq. (|3Tj) . Ci satisfies 



dr dr \r ipi 



In terms of x defined by r = itanx, it is rewritten as 



= — — (cos x + K sin x cos x + 3£/i sin x cos Xj • 



The integration of the third term is 



rx 

3i \i sin 2 x cos x^X = sin 3 X 

J-w/2 



sin x~r d Xi 
= ^/isin x+ / Co^sinx^X 

Jtt/2 

- ^sin 3 X - Co cos x 

where we have set ft, = 1 as explained in Appendix El Hence the solution is 

2k k 2 2k 

Ci = — (1- sinx) - -r cos 2 x+ -tCqcosx- 2/^sin 3 x. 



(35) 



(36) 



(37) 



(38) 



Since the radius of the photon sphere is r c — 1, X — r /£ an d the last term can be neglected. Hence we have 

Ci = (2(3 + 3k + 2)5, (39) 

where wc have introduced 5 = k/£. Moreover, from Eq. (|26|) . we find r c = 1 — 5. Putting these into Eq. (|25[) gives up 
to the first order 



a 2 + {f3 + 2k) 2 + 26(13 + 2k)(2/3 + 3k + 2) = 16(1 + 26). 



(40) 



Namely, the shadow is now the ellipse suppressed in the (3 direction and elongated in the a direction. If only one of 
a binary black hole is observed, the shape of its shadow determines not only its mass but also the information of the 
other black hole, 6 = k/£. The other black hole should exist in the direction where the shadow is suppressed. 




FIG. 7: The analytic solutions Eq. (|40|l (red) superimposed on the black hole shadow for l/m = 20, 10, 4 with m = ml . 



In Fig. [7j we show Eq. (|40| superimposed on the black hole shadows for several I with k = 1. Wc find excellent 
agreement. 
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FIG. 8: The black hole shadows for the two black hole system in the KT solution plotted in a-/3 space normalized by tM with 
each physical time at the observer tj\H\~ x = 0, 1.6, 3.2, 4.6, 7.6, 16.1. 



III. SHADOWS OF COLLIDING BLACK HOLES: KASTOR-TRASCHEN SOLUTION 

Finally, we study the shadows of colliding black holes in the Kastor-Traschen(KT) solution pfl. KT solution is 
a time-dependent generalization of the MP solution and describes an arbitrary number of extremely charged black 
holes in de Sitter universe. It is reduced to the Majumdar-Papapetrou(MP) solution when the positive cosmological 
constant A = and is reduced to the extremely charged (Q = M) Reissner-Nordstrom-de Sitter (RNdS) solution. 

The metric in the cosmological coordinate is given by 



ds = —a f2 dr + a £1 (dx + dy + dz ), 

1 / A „ N - rrii 

, H = ±\-, tt = l+y — 

Ht V 3 4^ an 



r,Ht 



(41) 
(42) 



where, r and t denote conformal time and physical time respectively. Here, H > {H < 0) corresponds to expansion 
(contraction). In the contracting universe (H < 0) the KT solution describes the collision of black holes. 

Let us consider the situation where an observer is near inside the cosmological horizon (r b s — > ?' + ) in the contracting 
coordinate. The cosmological horizon r + of Q = M RNdS solution in the cosmological coordinate is given by 



1 



2 H 



(1 + y/l-4M\H\) - M. 



We define the following parameters, which form the celestial coordinate system, as 



P(t) 



P(r) 



(43) 



(44) 



where PW are the momenta in the local inertial frame and ar b s is the physical distance between the observer and 
the center of the coordinate. 
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A. Shadows of Colliding Black Holes 



Let us consider a two black hole system as an example of colliding black holes. Each black hole is located at 
z = ±£/2 in the comoving coordinate. We set an observer at a fixed point inside a cosmological horizon in the 
physical coordinate. First we take 6 f, s = n/2 in terms of the polar coordinate. 

We then numerically calculate the null geodesies from the observer in the expanding coordinate. The null gccdcsics 
which eventually fall into the black hole horizons are regarded as shadows. Note that the time reverse is the null 
geodesies going from the black hole to the observer in the contracting universe. 

Fig. |S] shows that the shadows of two black holes with same masses m\ = m-i at each physical time t seen by 
observers at o b s = tt/2 with e = a\H\(r + — r Q b s ) = 0.01. We take M = m\ + ??i2 = 0.1/|i/|. The separation of two 
black holes is chosen as al = 4 x 10~ 3 /|if | at t = 0. Here, the celestial coordinates a and (3 are normalized by eM in 
order to keep the shape of the shadows independent of a location of the observer. 

At t = and t — 1.6/|Zf|, the black holes are mutually away enough. However, one can find that their shapes 
are a little bit elongated in the a direction and squeezed in the /3 direction from the circles with a radius of 
irriie/ yl + 4m; | H | ~ 1.82eM when they are considered as single black holes in a-/3 space. This deformation is 
caused by the existence of the other black hole in the opposite side as explained in Sec. II. D. 

At t = 3.2/\H\ (and even at t = 1.6/\H\), an eyebrow-like structure around each black hole appears. This kind 
of structure is quite unique to the multi black hole system. The reason why these structures appear is explained in 
Sec. II. C. If the impact parameter of the null geodesies is slightly smaller than the radius of the photon sphere, these 
geodesies will eventually fall into a black hole horizon. On the other hand, for a slightly larger impact parameter, 
the winding geodesies will gradually increase the distance to the black hole and eventually go away from the black 
hole, or fall into the horizon of the other black hole. The latter case creates the eyebrow-like shadow along the main 
shadow. The situation is quite similar to the orbits of the null geodesies in the Majumdar-Papapetrou solution. 




t=7.6 




Apparent horizon 



4 6 
t[x1/|H|] 
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12 



FIG. 9: The behavior of the null geodesies observed in the direction a = /? = 0ati = 7.6, 12.2, 16.1/|-ff|. The dashed curve 
is the evolution of the common apparent horizon enclosing the two black holes. The vertical axis is a(t)x. The null geodesies 
go through the middle of the black holes (a; = 0) at t ~ —2.2/\H\. Later, the black holes get closer and a common apparent 
horizon appears at t ~ 0.3/\H\. The null geodesies with t — 16.1 stays close to (but just outside) the horizon for a while and 
then reaches the observer when t — 16.1/1 -H" I . 



At t = 4.6/|iT| and at t = 7.6/\H\, the eyebrow-like structures grow and the main shadows come close each other. 
One can find there still remains a region where photons can go through between the main shadows. The reason why 
such a region remains is the following. In a single black hole system, a black hole horizon is enclosed with the photon 
sphere. On the other hand, in a two black hole system, two photon spheres intersect at the x-z plane where the null 
geodesies cannot fall into either one of the black holes. Accordingly the null geodesies can go through this plane, 
which corresponds to j3 = in the celestial coordinate until two black holes merge and form a horizon. Even at 
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t = 16.1/|-ff I, there still remains a region where photons can go through between the main shadows. 2 According to 
a common apparent horizon encompassing two black holes appears when at < 10~ 2 /|if|. Since we consider the 
null geodesies in the contracting universe, the geodesies coming into the observer at t can go near the black holes at 
the time much earlier than t when the distance between the black holes at is much larger than 0.01/|if| so that the 
geodesies can go through in between. In fact, some of the null geodesies observed at t = 16.1/|ff| pass through the 
middle of the black holes at t ~ —2.2/\H\ when the distance between the black holes is at ~ 0.036/|-ff|(> 10~ 2 /|-ff|) 
and there is no common apparent horizon (Sec Fig. [5]). 

Overall, the shapes of these shadows look quite similar to those in Fig. [2J 

B. Shadows of "Coalescing" Black Holes 

Finally, we consider the situation where one observes black holes from arbitrary azimuthal directions to mimic the 
coalescing binaries. We have calculated shadows for several different values of angle 6 b s at t = 3.7 in Fig. [10] As we 
decrease 9 b s from tt/2, the left main shadow of Fig. [5]becomes elongated, and eventually merges with the eyebrow-like 
structure of the right side and forms a ring structure surrounding the right main shadow. For comparison, in Fig. II H 
the shadows of the MP solution for several different angles o b s are shown. Again, we find that both look similar. 




-4 -2 2 4 -4 -2 2 4 

|3/(eM) 



FIG. 10: The shadows of two black holes in the KT solution at t — 3.7 /\H\ viewed by the observer at # b s = tt/2, n/3, tt/6, 0. 



IV. SUMMARY 

We have studied the null geodesies in the static/dynamic multi-black hole solutions: Majumdar-Papapetrou (MP) 
solution and the Kastor-Traschen (KT) solution. We have calculated the shadows of these multi-black holes and 
found that the shadows have structures distinct from the mere superposition of the shadow of each black hole: the 



2 The merger of the shadows found in is due to the low resolution of the numerical calculation. 
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FIG. 11: The shadows of two black holes in the MP solution with I = 2M viewed by the observer at # bs = t/2, 7t/3, 7t/6, 0. 

eyebrow-like structures outside the main shadows and the deformation of the shadows. We have presented analytic 
estimates of these structures using the MP solution to show that the width of the eyebrow is related with the distance 
between the black holes and that the shadows are deformed into ellipses due to the presence of the second black holes 
and the separation between the shadows is larger. 

These analytic results help us to have qualitative understanding of the features of the shadows of colliding black 
holes which are studied in our previous paper. We expect that following two features of black hole shadows are 
general and appear in more realistic situation. First one is the eyebrow-like structure which shows up during the 
merger process. Second is the deformation of the main shadow and the larger separation than the true distance. 

These features in the shadows can be used as probes to find the multi-black hole system at the final stage of its 
merger process. For that purpose, we have presented the shadows of the colliding black holes in the KT solution 
by changing the direction of the observer to mimic the coalescence of the binary black holes. In order to study the 
shadows of a realistic black hole binary, the effects of the accreation disk should also be considered, which is left for 
our future study. 
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Appendix A: Derivation of Eq. (1311) 

In this Appendix, we show that up to 0(1/ 1) ( satisfies Eq. (pTTl) . 
First we prove the following relation 



\p\<^ + 0{l/i). (Al) 
r 



For convenience, we introduce p = X 2 + £ 2 - First note that since C 2 — P 2 p 2 > from Eg . (f2"l"T) . we have p = for 
£ = 0. Therefore, it suffices to prove the relation for p ^ (or £ ^ 0). 
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Since p — > Q/r for r — > oo, from Eg. (1211) we have for r — > oo 



dp 
dr 



C / i-pV/C 2 



r 2 n 2 V l-V/(r 2 ft 4 )' 



(A2) 



The sign of the right-hand-side changes at p 2 p 2 = ( 2 or p 2 = r 2 Jl 4 . The latter corresponds to the turning point of 
the orbit (Eq. (fT9|0 . Before reaching the turning point, there can be several points r such that p 2 p 2 = ( 2 . Denoting 
the largest one among such r as r\ we first prove the relation for r > r\ . The relation is then consistent for r < r\ , 
because p takes a maximum value (/p at r\. 
The function ft defined by 



h = I i-pV/C 2 



1 -p 2 /{r 2 n 4 ) 



(A3) 



takes the minimum ft = at r = i\ and asymptotes ft = 1 for r — > oo (0 < ft < !)• Using ft, Eq. (|A2[) and Eq. 
become 



dp 
dr 
dc 
dr 



r 2 n 2 ' h ' 



Introducing g± = fi ± - , from Eq. (|A4[) and Eq. (|A5[) , g± satisfies 



The second term is 



n 

2—ip >fl h 



2£kQ 



Then the integration of the second term is 



(r 2 + I 2 - 2£rp) 3 / 2 



r n 



< A = - 2-i/vMr < - / 



ft >0. 



o 



(A4) 
(A5) 

(A6) 
(A7) 



-2k 
2k 

T 



tdr 



(r 2 + £ 2 ) 3 /2 



vV^+T 2 



1 + - + 0(1/^) 



C(l/£ 2 ), 



(A8) 



which shows that it is the first order quantity. 

Noting g± — > for r — > oo, the integration of Eq. (|A6[) gives 



(g + +A){g--A) 



Using Eq. (|A5p . the integrals in the right-hand-side of Eq. (|A9p can be rewritten as 



[/>" 


2 

+ 




-J oo ' 







(A9) 



r 1 



-%dr = ( / — dr + 2 / 0^,^^ 



—dr 



r c ft 



h 



dr + 2 Qipith 



di- 
ll 



-dr 



dr. 



(A10) 



The last term in Eq. (jAlOl) is less than A because of < ft/il 2 < 1, therefore this term is the first order quantity. 
Substituting these integrals into Eq. (|A9j) . we obtain 



(.9+ + A)(g_ - A) = -C 2 



r 




2 




'1 














7oo »" 2 . 









0(l/£). 



(All) 
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From < h/£l 2 < 1 we find that the first term in right-hand-side of Eq. (|A11|) is negative. Therefore, since 
A -0(1/*), 

g+g- < 0(1/*). (A12) 

Hence we have 

|M|<^ +0(1/4 (A13) 
r 

the relation holds at least up to the turning point of the orbit. 

Next, we show using the relation (|A1[) that Eq. (| A5|) can be simplified to give Eq. (pTTj) . Using the relation (|A1[) 
and Eq. (|A5[) . we have the inequality for r > p 



2 NhP) < § < - 2 ^W i-p 2 /(^ 2 ) - (AU) 



We denote r which satisfies r = p(r) as r p . Note that p = ^/ ( 2 + £ 2 ~ 0(1). Since the denominator in the square 
root at r = p is 1 — 1/S1 2 ~ 0(1), there exists constants c\ and C2 of 0(1), for < r p /r < 1 so that 



1+ ^>^V i-pvWr Vt^tt^)- 1 ^ 1 ^ 2 ' (A15) 



Then we show that the second term in the integral 



^(l + Cl (r p /r) 2 ), i-l,2, (A16) 

oo 



is of the second order. 

The second term becomes, neglecting higher order terms, 



(ra+ %/ a ( 1 + 1 /r)(r P /r) a dr 



2^ ^ sin X (l + cot X /t) cot 2 X dx ~2k-^ In ^, (A17) 



where we have introduced r = £tanx- Hence, the integral is of order 0(l/£ 2 ). 
For r < r p , from Eq. (|A5[) . we have 

C(r)-C(r p )^0(l/^ 2 ). (A18) 
Therefore, up to 0(1/*), Eq. (|A5[) is simplified by setting /i = 1 



dr 



This is Eq. ([3T 
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